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A quantum particle subjected to a constant force undergoes an accelerated motion follow¬ 
ing a parabolic path, which differs from the classical motion just because of wave packet 
spreading (quantum diffusion). However, when a periodic potential is added (such as 
in a crystal) the particle undergoes Bragg scattering and an oscillatory (rather than 
accelerated) motion is found, corresponding to the famous Bloch oscillations. Here we 
introduce an exactly-solvable quantum Hamiltonian model, corresponding to a general¬ 
ized Wannier-Stark Hamiltonian H, in which a quantum particle shows an intermediate 
dynamical behavior, namely an oscillatory motion superimposed to an accelerated one. 
Such a novel dynamical behavior is referred to as accelerated Bloch oscillations. Analyt¬ 
ical expressions of the spectrum, improper eigenfunctions and propagator of the gener¬ 
alized Wannier-Stark Hamiltonian H are derived. Finally, it is shown that acceleration 
and quantum diffusion in the generalized Wannier-Stark Hamiltonian are prevented for 
Airy wave packets, which undergo a periodic breathing dynamics that can be referred to 
as Airy-Bloch oscillations. 
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1. Introduction 


Stark Hamiltonians play a central role in several f ields o f physics, including quantum 
mechanics, condensed matter physics and optics i In non-relativistic quantum 
mechanics, the simplest case is that of a quantum particle subjected to a constant 
force F, which is described by the Hamiltonian 


^2 

H = -e-^ + Fx= epV + Fx 
ax‘‘ 


( 1 ) 


where the first term of the right hand side of Eq.(l), ep^, is the kinetic energy op¬ 
erator, the second term is the potential of the external force, and Px = ~'idx is the 
momentum operator. As is well known (see e.g. chap.8.1 of Ref.EI or Ref.Sl) the spec¬ 
trum of H is purely absolutely continuous with a set of (improper) eigenfunctions 
given by shifted Ary functions. Like in the classical limit, a quantum wave packet 
undergoes a uniformly accelerated motion but spreads owing to quantum diffusion. 
Indeed, after a change of reference frame (from th rest frame to an accelerated 
one) and a gauge transformation of the wave function the Stark Hamiltonian (1) 
is equivalent to the Hamiltonian of a freely moving particle, i.e. with F = 0 (see. 
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for example, I) . A much richer and physically relevant dynamics arises when a 
periodic potential is added to the external force, i.e. for the Hamiltonian 

H = epl + V{x) + Fx (2) 

where V (x) is the periodic potential with period d. Such a problem was originally 
studied by Bloch, Zener and Wannier in a series of seminal papers in connection 
with the mo tion of an electron in a crystal subjected to a uniform electric field 
| 8 | 9|10 | ll J12j assumption that a single lattice band is excited and the 

force is weak to neglect tunneling into other bands of the crystal (Zener tunneling), 
the quantum motion of the elect ron can be described by an effective single-band 
Wannier-Stark Hamiltonian 1211^. Remarkably, owing to Bragg scattering in the lat¬ 
tice a quantum wave packet in the crystal undergoes an oscillatory (rather than 
uniformly accelerated) motion with period Tg = 2TT/{Fd), the famous Bloch oscil¬ 
lations (BOs). While in natural crystals electronic BOs have never been observed 
owing to dephasing and many-body effects, BOs have been impressively demon¬ 
strated in a number of experiments after the advent of semicondu ctor su perlattices 
as terahertz radiation emitted from coherently oscillating electrons Quantum 

and classical analogues of BOs have been also proposed and observed in a wide 
variety of differen t physical systems, such as in ultracold ato ms and Bose- Ein- 
stein co ndensat es |16|17 | 18|19|^20 | 21|^ ^ optical waveguide arrays l^^l^4 | 25|26|27 | 28| 
photonic [^SEQ] and acoustic superlattices, metamaterials 1^, plasmonic systems 
|33|34 |^ quantum liquids to mention a few. There is a very extended lit¬ 
erature on BOs and their applications. Several authors h ave investigated BOs in 
presence of lattice disorder and inhomogeneities I36|37 | 38 |^ particle interacti on and 
correlation inhomogeneous fields SSl, comple x poten tials |47 | 48( 

time-dependent forces | 49|50| gQg gf non-classical states of light etc. As far 

as applications is concerned, BOs h ave been used as a tool to explore topological 
properties of band structures I 53 | 54 |^ j^ake p recisio n measurements of gravity 1^, 
and to realize beam splitters for matter waves ESiSTI 

In this work we introduce an exactly-solvable generalized Wannier-Stark Hamil¬ 
tonian in which the dynamics of a quantum wave packet shows a behavior interme¬ 
diate between the uniform acceleration of Hamiltonian (1) and the oscillatory BO 
dynamics of Hamiltonian (2). Namely, an oscillatory dynamics is superimposed to a 
parabolic motion, an effect that will be referred to as accelerated Bloch oscillations. 
The generalized Wannier-Stark Hamiltonian considered in this work has the form 

H = epl+ T{p^) + Fx, (3) 

where T{q) is a periodic function of q with period 27r/d and zero mean. It differs from 
Eq.(I) because of the additional term T{px) to the free-particle kinetic energy oper¬ 
ator. A possible physica l im plementation of such an Hamiltonian has been recently 
introduced in optics Ref.^^, where transverse light dynamics in a self-imaging opti¬ 
cal cavity is described by an Hamiltonian of the form described by Eq.(l). In Ref.^^, 
the limiting case e = 0 was mainly investigated, corresponding to pseudo-periodic 
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BOs dynamics and absence of quantum diffusion. Here we consider the more general 
case e 7 ^ 0 and derive analytical expressions of the spectrum, eigenfunctions and 
propagator of the Hamiltonian (3). In particular, we prove rigorously that for e 7 ^ 0 
the new phenomenon of accelerated Bloch oscillations arises: an oscillatory motion 
with period Tb = 27r/(Fd), i.e. analogous to BOs, is superimposed to the freely- 
accelerating motion of the wave packet (as if T = 0 in (3)). Such a result holds for 
any normalizable wave packet, but can be violated for an initial non-normalizable 
probability distribution. Interestingly, for an initial condition corresponding to a 
(non-normalizable) Airy wave packet, eigenfunction of (1), a periodic oscillatory 
dynamics without acceleration is found, a phenomenon that can be referred to as 
Airy-Bloch oscillations. 


2. Generalized Wannier-Stark Hamiltonian 

Let us consider a quantum system described by the one-dimensional Schrodinger 
equation 

i^=Hilj{x,t) (4) 

with Hamiltonian given by Eq.(3), where 

00 

n<i)= E Tn exp{indq) (5) 


is a periodic function of period /d with zero mean (Tq = 0) and Px = —idx is the 
momentum operator. Eor e = 0, the Hamiltonian (3) is a Wannier-Stark Hamiltonian 
with a periodic kinetic energy operator that is generally introduced in solid-state 
physics as an effective Hamil tonia n to describe single-band electron dynamics in 
slowly-varying external fields this case T{q) is the band dispersion curve, 

and the discretization x = nd {n = 0, ±1, ±2, ±3,...) at localized Wannier sites 
must be accomplished I2II31 As is well-known, for e = 0 and provided that the 
discretization x = nd is accomplished, H has a pure point spectrum given by an 
equally-spaced Wannier-Stark ladder with energy separation ujb = Fd = 2t:/Tb\ 
in real space an initially localized wave packet undergoes periodic BOs with period 
Tb- The generalized Wannier-Stark Hamiltonian (3) that we consider in this work 
differs fro m the effective Hamiltonian introduced in the theory of electron dynamics 
in crystals^ ^l^^l because (i) the variable x is continuous rather than discrete, and (ii) 
the term breaks the periodicity of the kinetic energy operator. In such a case, 
H has a pure absolutely continuous spectrum —oo<E<oo and the improper 
(non-normalizable) eigenfunctions 4>e{x)^ satisfying the eigenvalue equation 


H(j)E{x) = E(j)E{x) 


( 6 ) 
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can be analytically determined in terms of series of shifted Airy functions. The 
explicit expression of derived in Appendix A, reads 

(f)E{x)= ^ pnAi(f^^ {x-nd-E/F)j (7) 


where we have set 


oTZ f d 


Pn — 


2^ei/3^^i/6 


dq exp \ iqdn + ^ / . 


and the following normalization condition holds 

/ OO 

dx(j)*E,{x)4iE{x) = 6{E - E'). 

-OO 


( 8 ) 

(9) 


It is interesting to specialize the general result given by Eqs.(7) and (8) to the fol¬ 
lowing two limiting cases. 

(i) In the limit T{q) = 0, i.e. when the Hamiltonian (3) reduces to the standard 
Stark Hamiltonian (1), one has and thus the improper eigen¬ 

functions (7) reduce to the usual form (see®) 


(ii) In the limit e —>■ 0, i.e. when the Hamiltonian (3) has the form of the effective 
Wannier -Stark Hamiltonian of solid-state physics with a periodic kinetic energy 
operator taking into account that 


one obtains 


lim 

7—>-0 


Ai(a:/7) 

7 


S{x) 


OO 

4>e{x) = anS{x — nd — E/F) 

n— — oo 


( 11 ) 


( 12 ) 


where we have set 

^ n-K/d r ^ /‘Q 

= 2T,.pi/2 J ^^dqexp}^iqdn+ j d^T(^) j . (13) 

Equations (12) and (13) justify the result previously given in Ref.l^ and indicate 
that, even at e = 0, the spectrum of H is absolutely continuous. Such a result stems 
from the fact that the variable x here is considered a cont inu ous variable, rather 
than discretized as in single-band electron transport theory 11^. 
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3. Wave packet dynamics and accelerated Bloch oscillations 
3.1. Propagator 

The solution to the Schrodinger equation (4) with a given initial condition 
can be formally written as 

/ OO 

dyU{x,y,t)'ip{y,0). (14) 

-OO 

The propagator U, or its kernel U{x,y,t) entering in Eq.(14), can be determined 
from the spectral representation of H, i.e. one has 

/ OO 

dE(j)};{y)(j)E{x)exp{-iEt). (15) 

-OO 

Substitution of Eq.(7) into Eq.(15) and after some cumbersome calculations one 
obtains (see Appendix B for technical details) 


id{x,y,t) = (-iFxt ^ Gi{t)exp 

where we have set 

OO 

Gi{t)= ^ pnPn_iexp{iFtdn) 

n— — oc 


(^y — X + dl — 


Aiet 


( 16 ) 


(17) 


and pn are given by Eq.(8). Exact analytical expressions of the function Gi(t), and 
thus of the propagator (16), can be given in the important case of a sinusoidal 
shape of T{q), i.e. T{q) = KCOs{qd), which in electronic crystal theory describes 
a nearest-neighbor tight-binding lattice band. This important and exactly-solvable 
case is presented in Appendix C [see Eq.(C.6)]. 

It is interesting to specialize the form of the propagator given by Eq.(16) to the 
following two limiting cases. 

(i) In the limit T{q) = 0, i.e. when the Hamiltonian (3) reduces to the standard 
Stark Hamiltonian (1), one has Gi{t) = {e'^F)~^/^Spo and thus 


n{x,y,t) = 


1 f ^ 

exp I —ix xt — i —-— 


Airiet 


exp 


(y — X — eFt^y 
Aiet 


(18) 


which is precis^ the propagator of the Stark Hamiltonian (1) (see for example 
Eq.(7.4) in Ref.II^. 

(ii) In the limit e = 0, Eq.(16) shows a singular behavior and the propagator lA can 
be calculated as a limit of Eq.(16) as e —>■ 0. Since U is the kernel of an integral 
transformation [Eq.(14)], by use of the phase stationary method as e —0 one can 
set 


Tilexp (-F) 

V 47r*et \ Aiet J 


6{x) 


(19) 
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in all terms under the sum in Eq.(16). Using this result, after some straightforward 
calculations one obtains 

oo / oo \ 

U{x,y,t) = F exp{—iFxt) E a^_iexp{iFtdn)\ 6{x — y — dl), (20) 

1— — 00 \n — oo / 

where cr„ are defined by Eq.(13). 


3.2. Accelerated Bloch oscillations 

Propagation of an arbitrary initial wave packet tplXyO) is determined by the inte¬ 
gral transformation (14) with the propagator U defined by Eq.(16). To study the 
properties of wave packet evolution, let us distinguish three cases: 

(i) First case: parabolic motion. This case corresponds to the Stark Hamiltonian 



Fig. 1. Parabolic motion and spreading of a Gaussian wave packet for the Stark Hamiltonian (1). 
(a) Evolution of the probability density t)\^ on a pseudocolor map. (b) Evolution of the wave 
packet center of mass {x{t)) and width Ax. Parameter values are given in the text. 


(1), i.e. T{q) = 0. The propagator is given by Eq.(18), which describes particle 
acceleration and quantum diffusion. An initial wave packet undergoes uniform ac¬ 
celeration and the center of mass of the wave packet {x{t)) follows the parabolic 

trajectory ESI 

{x{t)) = {x{0)) + 2e{px{0))t - eFf. (21) 

Analytical expression for '4){x, t) can be given, for example, for an initial Gaussian 
wave packet distribution. In addition to the parabolic motion, wave packet spreading 
is observed as a result of quantum diffusion. An example of accelerated Gaussian 
wave packet is shown in Fig.I. In Fig.1(a) the evolution of the probability density 
\'ijj{x,t)\^ is depicted in a pseudo color map for parameter values e = 1/2, F = 0.2, 
r(g) = 0 and for the initial condition t/jiXjO) oc exp{—x^/w^) with w = 5. Figure 












time t 
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1(b) shows the evolution of the wave packet center of mass {x{t)) and width Ax(t) 
defined by 

/ OO 

dx x\^{x,t)\'^ (22) 

-OO 

f roc '1 1/2 

Ax{t) = Ij dx{x — {x{t)))'^\4’{x,t)\‘^ > . (23) 

Note that the wave packet center of mass follows the parabolic trajectory according 
to Eq.(21), and that wave packet spreading arises because of quantum diffusion. 

ill)Second case: pseudo-Bloch oscillations. This case corresponds to the limit e = 0 



space v space X timer 

Fig. 2. Pseudo-periodic oscillations of a Gaussian wave packet for the generalized Wannier-Stark 
Hamiltonian (3) with e = 0 and with a sinusoidal kinetic energy operator T{px) = K,cos(dpx). 
Evolution of (a) the probability density t)p, and (b) of the real part of the wave function 

Re(^(a:,t)) on a pseudocolor map. (c) Evolution of the wave packet center of mass {x(t)) and 
width Ax. Parameter values are given in the text. 


in Eq.(3) and displays a pseudo-periodic dynamics. In fact, at discretized times 
t = tjn = utiTb, integer multiplies than the BO period Tb = 27r/(Fd), from Eq.(20) 
one obtains 

n{x,y,tjn) = exp{-iFxtm)6ix - y) (24) 

where we used the property crn<Xn-i — Hence from Eqs.(14) and (24) 

on has 

ipix^tm) = e'xp{—2'Kimx/d)'ip(x,Q). (25) 

Equation (25) shows that the dynamics is pseudo-periodic, i.e. the probability dis¬ 
tribution |'i/;(a:,t)p undergoes a periodic dynamics with the BO period Tg, i.e. 
\tp{x,tm)\^ = |■^/'(a:,0)p, whereas the wave function ilj{x,t) does not [this is due 
to the phase factor exp(—27r/ma:/d) in Eq.(25)]. An example of pseudo-periodic dy¬ 
namics for an initial Gaussian wave packet is shown in Fig.2 for a sinusoidal shape 
of T{q), i.e. T{q) = KCOs{qd). Parameter values used in the simulations are e = 0, 
F = 0.2, K = 1 and d = 4, corresponding to a BO period Tb = 27r/(Fd) ~ 7.85. The 
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initial condition is the Gaussian wave packet 0) oc exp(—with w = 5 as 
in Fig.l. Figure 2(a) and (b) depict on a pseudo color map the temporal evolution 
of the probability density |'0(a:,t)p and of the real part Re('0(x,t)) of the wave 
function, respectively. Note that, while the probability density undergoes a periodic 
oscillatory behavior analogous to the ordinary BOs in a crystal, the wave function 
does not; in particular the interference fringes visible in Fig.2(b) arise from the 
phase factor exp(—jFxt) appearing in Eg.(20). Such a pseudo-periodic dynamical 
regime was previously introduced in Ref.'^ and referred to as pseudo-Bloch oscil¬ 
lations. 


(iii) Third-case: accelerated Bloch oscillations. This is the most general case, corre¬ 
sponding to T{q) 7 ^ 0 and e 7 ^ 0. In this case the dynamical behavior is intermediate 
between the two previously considered regimes, i.e. one observes an oscillatory dy¬ 
namics of the wave packet superimposed to the parabolic path as described by 
Eq.(2I). Such a property can be proven by observing that, at times tm = itiTb 
(m = 0, ±1, ±2,...), the two propagators W(a:, y, t) as given by Eq.(I 6 ) and Eq.(18) 
do coincide. In fact, at times t = tm one has 

CXD - 

Gl{tm) = PnPn-l = = 2/3pl/3^*'0 


where fl; are defined by Eq.(A.I7) in Appendix A. Substitution of Eq.(26) into 
Eq.(I 6 ) yields 


Gix,y,t) = 


1 


Airietn 


exp ( —iFxtm — i- 


eFH: 


2^3 


exp 


jy-x-cFtl^Y 

^ietm 


(27) 


which coincides with Eq.(I8) taken at t = tm. Indicating by Fli and 772 the two 
Hamiltonians defined by Eqs.(l) and (3), respectively, i.e. Hi is the limiting case of 
H 2 when T{q) = 0, the previous result can be formally written as^^ 

exp 772 ) = exp(-i<„77i) (28) 


i.e. if the dynamics is mapped at discretized times tm integer multiplies than the BO 
period Tb the two Hamiltonians (I) and (3) yields the same evolution. Therefore 
the oscillatory motion in each BO cycle is superimposed to the parabolic path (21). 
Such a dynamical behavior can be thus referred to as accelerated Bloch oscillations. 
At times different than tm = uiTb, the following relation between the solutions 
ij}i{x,t) = exp{—itHi)'ip{x,0) and 'tp 2 {x,t) = exp{—itH2)'ijj{x,0) to the Schrodinger 
equations with Hamiltonians Hi and H 2 corresponding to the same initial condition 
ip{x,t) can be derived (see Appendix D) 

00 

il) 2 ix,t) = ^ Ai{t)tpi{x - ld,t), 

I —— 00 


(29) 
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where we have set 

OO 

Ai{t) = F^^^e^^^Gi{t)exp{—iFdlt)=F^^^e^^^ ^ /o*/9„+; exp(iFdtn). (30) 

n— — oo 


For example, for a sinusoidal function T(q) = KCOs{qd) the explicit expressions of 
the functions Ai{t), as obtained from Eqs.(30) and (C.5), read 

A,(t) = J.d4,V 2-2co.{Fdt) ) 


\Fd 


1 — exp{iFdt) 


exp(—zFdZt). (31) 


Equation (29) indicates that, at times t ^ tm, the solution 'ild 2 {x,t) is given by a 
suitable superposition (interference) of shifted replica of ipi(x,t), weighted by the 
complex amplitudes A/(t). At t = tm, Ai{tm) = and thus ip 2 {x,tm) = f’lix, tm) 
as previously discussed. 



space X t 

Fig. 3. Accelerated BOs of a Gaussian wave packet for the generalized Wannier-Stark Hamiltonian 
(3) with a sinusoidal kinetic energy operator T{px) = K.cos{dpx). (a) Evolution of the probability 
density \'ilj{x,t)\^ in a pseudo color map, and (b) of the wave packet center of mass and width. 
Parameter values are given in the text. 


An example of accelerated BOs is shown in Fig.3 for a sinusoidal function T{q) = 
KCOs{qd) and for parameter values e = 1/2, F = 0.2, k = 1 and d = 4, i.e. for the 
same parameter values as in Fig.2 except for e = 1/2. The initial condition is the 
Gaussian wave packet ip{x,Q) oc exp(—x^/rc^) with w = 5 (i.e. as in Figs.l and 
2). The figure clearly shows that the wave packet undergoes and oscillatory motion 
with period Tb over the averaged parabolic path of Fig.l. Note that, while in the 
pseudo BOs regime [case (ii), see Fig.2] wave packet spreading is suppressed, in the 
regime of accelerated BOs (Fig.3) the wave packet spreads on average following the 
same spreading law of the parabolic case (i) discussed above. 

4. Airy-Bloch oscillations 

As shown in Sec. 3.2, any normalizable wave packet undergoes accelerated BOs and 
spreading when evolved by the generalized Wannier-Stark Hamiltonian (3) with a 
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non-vanishing value of e. However, such a property can be violated by an initial 


wave packet that is not normalizable, i.e. such that f^^dx = oo. An 


interesting case is the one corresponding to an initial condition tjj(x,0) which is a 
generalized eigenfunction of the Stark Hamiltonian (1), i.e. an Airy wave packet. 
Let us assume, for example, the eigenfunction of (1) with energy E = 0, which apart 
from a normalization constant is given by [see Eq.(lO)] 



(32) 


For the free-particle Schrodinger equation, i.e. for the Hamiltonian H = ep^, Airy 
wave packets are shape preserving ones and un dergo a self-accelerating motion, as 
originally shown by Ber ry a nd Balazs in Ref.^^ and extended in several subsequent 
works (see, for example, and references therein). For the Stark Hamiltonian (1), 
they do not accelerate, i.e. they are at rest, because of the additional force F: indeed 
they are eigenstates of the Hamiltonian (1). Interestingly, we show now that for 
the generalized Wannier-Stark Hamiltonian (3) the Airy wave packet (32) evolves 
undergoing a periodic breathing dynamics with the BO period Tb = 2TT/(Fd). Such 
a periodic and acceleration-free breathing dynamics is referred to as Airy-Bloch 
oscillations. In fact, the evolution of the wave packet (32) under the generalized 
Wannier-Stark Hamiltonian (3) can be readily calculated using Eq.(29) and taking 
into account that '0(x, 0) is an eigenfunction of (1) with zero energy. One obtains 



(33) 


Since for tm = ttiTb, one has ip{x,tm) = ipiXjO), i.e. a strict peri¬ 

odic dynamics is obtained. Hence for the initial Airy distribution (32) of the wave 
function the acceleration motion is suppressed, and a periodic breathing dynam¬ 
ics is established. An example of Airy-Bloch oscillations is shown in Fig.4(a). The 
periodic dynamics shown in the figure resembles a quantum carpet^^ in Hamilto¬ 
nians with absolutely continuous spectrum, such as quant um car pets arising from 
the Talbot effect for the free-particle Hamiltonian H = ep^ | 63 | 64 | 

In practice, the Airy distribution is an idealized one and corresponds to a delo¬ 
calized (non-normalizable) state of the particle. However, it can be approximated 
by normalizable distributions obtained by enveloping the Airy function with a 
sufficiently-decaying function at x —>■ —oo. For example, an initial normalizable 
wave packet distribution that approximates Eq.(32) is given by 



(34) 


where a > 0 and V is a normalization constant. Note that the function defined 
by Eq.(34) goes to zero at both x —>■ ±oo and is normalizable because the Airy 
function decays at x —>■ oo faster than exponential. Since the wave packet defined 










time t time t 
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space X 


space X 


time t 


Fig. 4. Airy-Bloch oscillations, (a) Quantum carpet, showing the evolution of the probability 
density \'ip{x,t)\^ in a pseudo color map for the initial condition (32) (ideal Airy wave packet). 
Parameter values are e = 1/2, F = 0.2, T{q) = K.cos{qd) with d = A and k = 1. (b-c): same as 
(a), but for the initial wave packet distribution (34) and for increasing values of a: a = 300 in (b), 
a = 100 in (c) and a = 50 in (d). (e) Evolution of the revival probability [Eq.(35)]. Curves 1,2 and 
3 refer to the simulations in panels (b), (c) and (d), respectively. 


by Eq.(34) is normalizable, its center of mass undergoes an oscillatory and (on 
average) accelerated motion; however, for a sufficiently large values of a the Airy- 
Bloch breathing dynamics of Eig.4(a) can be observed in the earlier BO cycles, as 
shown in Fig.4(b-d). After a few BO cycles, the periodicity is lost according to the 
analysis of Sec.3.2. This can be seen by computing the temporal evolution of the 
revival probability into the original state, defined as 


rrev(i) — 


dx 'tp*{x, 0)ip{x, t) 


(35) 


The behavior of Prev{t) is shown in Fg.4(e). Note that, as a decreases, periodicity 
of the dynamics is rapidly lost. 


5. Conclusion 

In this work the phenomena of accelerated and Airy-Bloch oscillations have been 
predicted, which provide significant extensions of the famous Bloch oscillations 
originally predicted for electrons in a crystal under a uniform electric field. We 
introduced an exactly-solvable generalized Wannier-Stark Hamiltonian and showed 
rather generally that a dynamical regime intermediate between a pure oscillatory 
and an uniformly accelerated motion arises (accelerated Bloch oscillations). As a 
special case, for wave packets with an Airy shape acceleration can be suppressed 
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and a pure oscillatory (breathing) dynamics, leading to quantum carpets, is pre¬ 
dicted (Airy-Bloch oscillations). Owing to the possibility to emulate Schrodinger 
equat ions with engineered potentials and kinetic energy operators offered by optics 
| 65 | 66 |^ predicted phenomena of accelerated and Airy-Bloch oscillations could be 
observed in an optical setting, as discussed in Ref.lSS]. 


Appendix A. Generalized eigenfunctions of the Hamiltonian (3) 

In this Appendix we derive Eqs.(7) and (8), given in the text, for the improper 
eigenfunctions of the Hamiltonian (3). To this aim, let us first note that for T^) = 0 
the eigenfunctions of H are simply given by shifted Airy functions (see e.g. For 
T{q) 7 ^ 0, since T{q) is periodic the action of the operator T{px) on a generic 
function f{x) is to sum up shifted replica of the function, namely one has 

T{px)f{x) e:xp{nddx)f{x) = '^Tr,f{x + nd). (A.I) 

n n 

Hence we may look for a solution to the eigenvalue equation (6) of the form 

OO 

(j)Eix) = ^ pnAi{a{x - nd-13)), (A.2) 


where the coefficients p„, a and (3 in Eq.(A.2) are unknown parameters at this stage. 
Substitution of the Ansatz (A.2) into Eq.(6) and taking into account Eq.(A.l) and 
that d^k\{x)/dx^ = a:Ai(x), one obtains 

{E — Fx) PnAi (a{x — nd — f3)) = —ea^ — nd — /3)pnAi [a(x — nd— /3)) 

n n 

+ '^^TipnAi{a{x — nd + Id — 13)). (A.3) 

1,71 


To satisfy Eq.(A.3), the coefficients of the terms multiplying x on the left and right 
hand sides of the equation must coincide. This yields the following expression for 
the coefficient a 



and Eq.(A.3) simplifies as follows 


(A.4) 


E PnAi {a{x — nd — (3)) = F'^^{nd + j3)pnA\ {a{x — nd — j3)) 

n n 

+ Tipn+^ Ai (Q!(a: - nd - (3)) . (A.5) 


Equation (A.5) is satisfied provided that the following difference equation among 
the coefficients /?„ holds 

{E - l3F)pn = Fdnpn + ^ Tipn+i- 

i 


(A.6) 
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To solve Eq.(A.6), let us introduce the Fourier spectrum 

OO 

S{q) ='^ Pn ex.p(-iqdn) (A.7) 

— CXD 

which is a periodic function of q with period 2Trjd. The coefficients are determined 
from he spectrum S{q) by the inverse relation 

d 

Pn = 7r dqS{q)e^p{iqnd). (A.8) 

J_7r/d 

From Eqs.(A.6) and (A.7) the following differential equation for the spectrum is 
obtained 


^F— = [E-|3F-T{q)]S{q) (A.9) 

which can be solved, yielding 

S{q) = S{0) exp |^(F; - pF)q + ^ £ d^O^ ■ (A.IO) 

The parameter (3 is determined by imposing the periodicity of S{q), i.e. that S{q + 
27r/(i) = S{q). Since dqT{q) = Tg = 0, the periodicity condition on S{q) is 
ensured by assuming 


P = E/F (A.ll) 

so that one has 

S{q) = SiO)e^pi^^ (A.12) 

Substitution of Eq.(A.12) into Eq.(A.8) yields 

Pn= ^ 

Finally, the constant S'(O) is determined by imposing that the generalized eigenfunc¬ 
tions (j)E{x) satisfy the usual normalization condition (9) given in the text. To this 
aim, let us explicitly calculate the scalar product {(j)E'{x)\(j)E{x)) using Eq.(A.2). 
One obtains 

{4)e’{x)\4)e{x)) = p*^pi{A\{a{x - nd - E' / F)\k\{x -Id - E/F)). (A.14) 

n,l 


Taking into account that 

/■°° 1 

/ dxAi[a{x — a)]Ai[Q!(a; — &)] = —d(a — b) 

J —OO ^ 


(A.15) 
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one has 

{(I>e'{x)\(1)e{x)) = -^ PnPi^i^' - E + Fdn - Fdl) = 

n,l 

= ^y"niS(E'-E-Fdl). 

i 


where we have set 


(A.16) 


ni=Y,P*nPu+l (a.17) 

n 

Using Eq.(A.7) it can be readily shown that 

ni =-^ [ dq\S{q)\‘^ exp{iql). (A.18) 

^ 7 ’" J-n/d 

Since |5'(g)p = |5'(0)p [see Eq.(A.12)], one has U; = |S'(0)p(I;^O; and thus 

{(^E'ixMEix)) = _ E). (A.19) 

Hence the normalization (9) is obtained by assuming 

cx 1 

>5'(0) = (A.20) 

Substitution of Eq.(A.20) into Eq.(A.13) yields the expression (8) given in the text 
for the coefficients p„. 


Appendix B. Propagator of the Hamiltonian (3) 

In this Appendix we derive the analytical form (16) given in the text for the kernel 
U{x, y, t) of the propagator tj{t) of the Hamiltonian (3). Substitution of Eq.(7) into 
Eq.(I5) yields 


U{x,y,t) = '^^Pl Pn / dEAi{a{y — Id — E/F)) Ai {a{x — nd — E/F)) ex.p{—iEt) 

n,i 

(B.l) 

with a = (E/e)^/^. After a change of the integration variable and summation 
indices on the right hand side of Eq.(B.l) one obtains 


U(x, y,t) = — exp{—iEtx) PnPn-i exp{iFtdn) i $(?/ — x + dl, t) (B.2) 

i \ n J 

t 

/ OO 

(i^Ai(^)Ai(^ + ax) exp{iFt^/a). (B.3) 

-OO 


where we have set 


The integral of the product of Airy functions on the right hand side of Eq.(B.3) 
with a given by Eq.(A.4) can be calculated in a close form and reads 




. / 1 

(x - eFt^) ieF'^t^ 

V dTTzeto^ 

Aid. 3 


(B.4) 
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After setting 


one then obtains 


Gi{t) = ^ PnPn-i exp{iFtdn) 


(B.5) 




2 ,3 \ °° 


exp I —iFxt — i 




^ Gi(t)exp 
1 — — 00 


{y — X + dl — eFt^^ 


Met 


(B.6) 


From Eqs.(A.4) and (B.6) then follows Eq.(16) given in the text. 


Appendix C. Propagator for a sinusoidal lattice band 

Let us consider the case of a sinusoidal shape for the function T(q), namely 

T{q) = KCOs{qd) (C.l) 

which corresponds in ordinary theory of electronic states in crystals to the well- 
known tight-binding lattice model in the nearest-neighbor approximation. In this 
case from Eq.(8) one obtains 

' - (-—) . 

V Fd) 


In deriving Eq.(C.2) the following integral representation of Bessel J„ function has 
been used 


1 r 

Jn{y) = — / dx exp{inx — iy sin x). 
27 ^ J-TV 


(C.3) 


The functions Gi{y), given by Eq.(17), take the form 
1 "" 


= ,2/3fi/3 E Jn{-^)jn-i(-^)exp{iFtdn). (C.4) 


n— — cio 


The sum on the right hand side of Eq.(C.4) can be calculated using the Graf’s 
addition formula for Bessel functions ( cf. Sec. 11.3 (1) of Ref.yielding 


= g2/3fi/3 -^-^ (^^^/2-2cos{Fdt)) 


1 — exp{—iFdt) 
1 — exp{iFdt) 


- 1/2 


(C.5) 


Finally, the propagator (16) takes the analytical form 

exp ( —iFxt — ^-^^2 — 2cos(Fdt)^ x 


U{x,y,t) = 


Mxiet 


1 — exp{—iFdt) 

1/2 

— X — dl — eFt^^ 

1 — exp(iFdt) 


Met 


(C.6) 
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6. Proof of Eq.(29) 

Let us indicate by 

, ■i/' 2 (a:, t) = exp(—0) (C.7) 

the solutions to the Schrodinger equation idt'^p = H%l) with Hamiltonians Hi = 
ep^ + Fx and H 2 = Hi + T{px), respectively, corresponding to the same initial 
condition '4){x, 0). One has 

POO 

'ipi(x,t)= dyUi{x,y,t)'ip{y,0) (C.8) 


ip 2 {x,t)= dyU2{x,y,t)'tp{y,0) (C.9) 

J —00 

where the propagators Ui{x,y,t) and U 2 {x,y,t) are given by Eqs.(18) and (16), 
respectively. Erom a comparison of Eqs.(16) and (18) one can readily write 

00 

U 2 {x,y,t)= ^ Ai{t)Ui{x - ld,y,t) (C.IO) 

I —— CO 

where the functions Ai{t) are defined by Eq.(30) given in the text. Substitution of 
Eq.(D.4) into Eq.(D.3) finally yields Eq.(29) given in the text. 
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